Boundary Value Problems in MATLAB